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Abstract. Polycrystalline thin films can be unstable with respect to island formation 
(agglomeration) through grooving where grain boundaries intersect the free surface 
and/or thin film-substrate interface. We develop a phase-field model to study the 
evolution of the phases, composition, microstructure and morphology of such thin 
films. The phase-field model is quite general, describing compounds and solid solution 
alloys with sufficient freedom to choose solubilities, grain boundary and interface 
energies, and heats of segregation to all interfaces. We present analytical results which 
describe the interface profiles, with and without segregation, and confirm them using 
numerical simulations. We demonstrate that the present model accurately reproduces 
the theoretical grain boundary groove angles both at and far from equilibrium. As 
an example, we apply the phase-field model to the special case of a Ni(Pt)Si (Ni/Pt 
silicide) thin film on an initially flat silicon substrate. 



PACS numbers: 81.10.Aj, 07.05.Tp 

Keywords: phase-field method, polycrystal, grooving, segregation 

Submitted to: Modelling Simulation Mater. Sci. Eng. 
1. Introduction 

While grain boundary grooving as a source for agglomeration in thin films is well known 
PP, most theoretical models for grooving- induced agglomeration are too idealized to 
be useful to predict agglomeration in technologically interesting materials. Most such 
materials are alloys containing many elements and often multiple phases. In such 
materials, segregation of one or more of the alloy components to one or more of the 
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interfaces (surface, grain boundaries, film/substrate interface) is the rule, rather than 
the exception. Hence, models based upon grain boundary grooving in single component 
films with no intermixing between the film and substrate are useful to understand the 
phenomenon of grain boundary grooving-induced agglomeration but cannot be relied 
upon to predict the rate of grooving or the phenomenology that has its origin in 
alloying effects. Unfortunately, including the full richness of such phenomena in real 
materials necessarily forces us to employ numerical, rather than analytical, treatments. 
In this paper, we describe a phase-field-based method for modelling the evolution of 
grain boundary grooves that is sufficiently general to include interfacial segregation, 
compound formation, and interdiffusion. 

The grain boundary grooving phenomenon has been known and studied 
theoretically for some time [21 El IH IHl IE]- Grooving is a general phenomenon which 
occurs at the free surface of all polycrystalline materials. While the model for grain 
boundary grooving described below is quite general, for specificity we focus on a 
particular example: a polycrystalline metal silicide film on a silicon substrate. The 
formation and evolution of thin polycrystalline metal silicides or germanides is key to 
the performance of many semiconductor devices. However, these films can be unstable 
at high temperature through agglomeration into silicide islands via grain boundary 
grooving 7]. The films may then lose connectivity. Experimental observations suggest 
that alloying NiSi with Pt [HI Ej or implanting BF2 [HI HHj delays or suppresses NiSi 
agglomeration. Since segregation of platinum or fluorine may play a role in this 
phenomenon, our model will include the possibility for one or more species to segregate 
to a grain boundary and/or interface. 

The phase-field approach has been extensively used for predicting microstructure 
evolution in different materials processes; for example, solidification, ferroelectric 
transformations in thin films, spinodal decomposition, grain growth, and dislocation 
dynamics IT^ Ej- Since our focus in the grain boundary grooving/thin film 
agglomeration modelling is on a scale that is large compared to the interatomic distance, 
yet small compared with the grain size, the phase-field method provides an excellent 
means through which to model the microstructure/morphology evolution. 

In this paper, we develop a phase-field model to describe the evolution of the 
microstructure and morphology of a polycrystalline thin film alloy, including atomic 
diffusion, phase transformation, grain boundary grooving and morphology evolution. 
To this end, we construct a general free energy function that allows for an arbitrary 
number of grains, phases and atomic constituents. Analytical results for the profiles 
of the order parameters at a grain boundary are presented within the framework of 
this free energy model. We also discuss the numerical procedure used for evolving the 
phase fields. Finally, we present a series of results on grain boundary grooving and Pt 
segregation in a polycrystalline Nii_a;Pta;Si film on Si. 
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Figure 1. A scliematic illustration of a polycrystalline NiSi film on a silicon substrate. 
The are the order parameters. The definition of the groove angle 6 is also shown. 



2. Description of the phase-field method 

We consider the case of a polycrystalhne film above which is a vacuum and below which 
is a substrate. Figure ^ is an illustration for the case of a NiSi film on Si. The model 
presented below also allows for the formation of other phases (e.g., NiSi2). 

To describe grain boundary grooving and atom diffusion during microstructure 
evolution (we also allow for phase transitions), we use two sets of phase-field variables 
in our model. One set is order parameter fields which describe phases, while the other set 
describes composition fields. In order to describe compounds, we distinguish between 
different sublattices. We define c'p as the concentration of species j on sublattice i. 
In the case of the Ni-Pt-Si-Ge system, sublattice 1 refers to the metals, which may 
substitute for one another, and sublattice 2 refers to the group-IV elements: on the 
metal sublattice c^P and c^^^ are the concentrations of Ni and Pt and on the group-IV 
sublattice Cg^"* and C2^^ are the concentrations of Si and Ge. The order parameters ifjk 
indicate the fraction of phase k present. Each phase k may be polycrystalline, 
represents grain m of phase k. In what follows the index i will always refer to types of 
atoms (i.e., sublattices), j to species, k to phases, and m to grains. 

In phase-field theory, the rate of evolution of the system is related to the variation 
of the total free energy of the system. The evolution of the conserved fields (i.e., 
compositions) is routinely described through the Cahn-Hillard equation jT3j 



^ = VMV- 



(1) 



dt 5c{r,t)' 
The non-conserved order parameter fields are governed by the Allen-Cahn equation [T 
di,{v,t) _ 5G 

dt S^{r,t)' ^ ' 

In these equations G is the total free energy of the system, c is a concentration (conserved 
field), t/) is a non-conserved field (phase), and M and L are atom and interface mobilities 
respectively. 5G/5c and SG/Sip are functional derivatives of G with respect to c and 
ip, respectively. 
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3. Total free energy 



In order to apply (P) and (j21), we first need an expression for the free energy of the 
system that includes all of the physics necessary to describe the agglomeration process 
in polycrystalline thin films. Such free energies have been obtained for many different 
physical situations [ini HZl UHl HH] • The present derivation of the free energy is based 
upon some of these earlier constructions. 

Ignoring elastic contributions to the free energy for now, the total free energy of 
the system can be expressed as 



G'(W,{c})= / MW,{c})-s(W,{c})r + 



V 



(i) 



Va 



(i) 



dV. 



(3) 



In this expression, T is the temperature and h{{ip},{c}) and s ({■?/'}, {c}) are the 
enthalpy and entropy of the system, respectively. The integral is over the entire system. 
While all grains of the same phase k exhibit the same thermodynamic behaviour, each 
type of atom is unique. This explains why k depends only upon the phase k, while A 
depends on both the sublattice i and the atomic species j. 

We first focus upon uniform phases (i.e., the gradient terms in are set to zero). 
In each bulk phase, one of the order parameters is unity, while those describing all other 
phases are zero. To ensure this, we include a term in the enthalpy, hi, that has minima 
at ipk"'' = and 1 for each k and m. As is commonly done, we use a double-well 
potential to describe this contribution to the enthalpy: 



It is also necessary to ensure that the phase of interest corresponds to a minimum 
in enthalpy at a particular stoichiometry. In our reference application, we consider four 
phases. Si, MSi, MSi2 and vacuum, where M represents the metal species (either Ni or 
Pt or some combination of these). In the MSi2 phase, a third of the sites are occupied by 
metal atoms: CNi + cpt = 1/3. Therefore, the enthalpy of this phase must be a minimum 



at 



Hj(A;) = 1/3 for i = metals and k = MSi2. We describe the enthalpy as a 



quadratic function for each type of species with a minimum at the desired stoichiometry: 



(m) 




where Ci is the stoichiometry of the system at any point in space: 



(5) 



(6) 



^mV^i™'' is amount of phase k and Ei{k) represents the stoichiometry of phase k. 
The second parenthesis in (0) is the difference between the local concentration of atoms 
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Figure 2. The enthalpy of a Ni-Si system with two phases, (Si) and NiSi, as a function 
of CNi and V'NiSi- 



of type i and stoichiometry. Ai{k) is the solubility of atoms of type i in phase k. Large 
Ai{k) implies a small solubility of i atoms on the "wrong" sublattices in phase k. 

Although the double-well potential of (j^ ensures that each is either or 1 in 
the bulk phase, it is possible for the phase fractions {ip) of several different phases to be 
1 simultaneously (e.g., simultaneously 100% monosilicide and 100% vacuum). Clearly, 
this is unphysical. To rectify this problem, we add a contribution to the enthalpy that 
contains cross terms between the order parameters: 



where the coefficients X^y are chosen to be sufficiently large so that the total phase 
fraction at any point is close to unity. A similar expression is used for grain boundaries 
(i.e., where k = k'): 



Figure 121 shows that the enthalpy of a Ni-Si system with two phases — (Si) and NiSi — 
has a minimum corresponding to the silicon phase at (cni = 0,^NiSi = 0) and another 
corresponding to the NiSi phase at (cNi = l/2,'?/'NiSi = 1). 

There is a contribution to the free energy associated with the entropy of mixing of 
the different atomic species on each sublattice. The configurational entropy Si associated 
with occupancy of sublattice i is 








Phase-field model for GB grooving 



where ri'p is the number of atoms of species j on sublattice i. We can rewrite ri'p in 



terms of the concentrations cP^ as n, 



Si) 



Ni/Ci, where Ni is the total number of 
atoms of type i. We can then rewrite this as the entropy per unit volume is the 
atomic volume on sublattice i) and expand it using the Stirling approximation, as 



k 



Si) 



Si) 



^ V In 



(9) 



The total free energy of the system is, therefore. 



G ({#■}, {c}) 



^ k m i i 



Vd 



(i) 



+ 



(m) 
k 



E4 



a. 



5Z 5Z 



k k'<k 




EE 

. m m'<m 



+ 



^B^E^E^'^'''^- 



+ 



(10) 



An elastic energy term can also be included in the free energy in order to account for 
possible elastic misfit between the various phases as shown in pOj l2T]. 

Interfacial segregation is a common feature in most multiphase and/or 
polycrystalline materials. Segregation of an element to the interface occurs if the free 
energy of the system is lowered by moving that element from the interior of a grain or 
single phase region to the interface. In order to incorporate such an effect, we modify one 
of the contributions to the free energy expression that controls interface energies. 
The appropriate enthalpy contribution is (interfaces) or (grain boundaries). This 
is accomplished by replacing X^.^/ with Xf^^i — Ykk'{i,j) c[''\ The set of parameters 
Ykk'ihj) describes the lowering of the k k' interface (or the grain boundaries of phase k 
if /c = k') energy upon segregation of species j (on sublattice i). 



4. Analytical results 

Before presenting a series of phase-field simulation results for the evolution of the 
polycrystalline silicide thin film, we examine how the general model described above 
describes segregation in order to confirm that the model is capable of reproducing the 
correct physics. In particular, we examine the case of a grain boundary between two 
semi-infinite grains, ^(1) and ^(2) 

in one dimension. Since only two grains of the same 
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Figure 3. -0 (left) and c (right) as a function of x/W . The hnes are obtained from 
(|12|l and (|14ll and the circles numerically. 



phase are considered the only unknowns are the phase fraction of grain 1, ip{x), and the 
composition profile across the grain boundary. We assume that there is only one species 
per sublattice except on one sublattice where there are two elements, one of which 
segregates to the grain boundary. As an example, we consider the case of Ni(Pt)Si. 
Hence, all compositions are fixed in terms of the Pt composition profile, c{x). The total 
energy is 

G'(V',c)= / (2D + XGB-VGBc)^2(l-^)2 + (ANi + Ap0c'2 + 

Jv 

2k i)'^ + A;bT [(1 - 2c) ln(l - 2c) + 2c In 2c] dV, (11) 

where the primes indicate derivatives with respect to the space coordinate x. 

If the segregation parameter Yqe is zero, the free energy can be split into two 
integrals, one depends on ip only and the other on c only. Therefore, equilibrium can 
be found by minimizing the contributions from ip and c independently. The solution to 
5G/5^ = is [131 

^(x) = i(l + tanh|p), (12) 
where W is the width of the grain boundary: 

= — (13) 
2D + Xgb ^ ^ 

This is the canonical profile of an interface where there is no interfacial segregation. 

Since Ygb = 0, the concentration profile c{x) is uniform. 

If Ygb 7^ 0, the concentration profile becomes 

" + 2l 2C + xT- Fob c„ A„, I Ap. W 
where Cq is the nominal concentration. The derivation (and assumptions) is presented 
in the appendix. Figure El shows ipi^) and c{x) obtained from (fT^ and (fT^ and from 
a direct numerical solution of the coupled set of equations. The results indicate that 
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grain boundary segregation does occur, its width is determined by the grain boundary 
width W. Figure El shows very good agreement between p4|) and simulations. 

5. Phase-field simulation results 
5.1. Numerical method 

In this section we show simulation results for two-dimensional systems with periodic 
boundaries along both directions. The evolution equations, (|T)) and 0, are solved in 
reciprocal space using a semi- implicit method [22] • 

The direct implementation of the phase-field method is complicated by the fact 
that the compositions and order parameters are not a priori to lie within the physically 
meaningful range, i.e. between zero and unity. To address this issue, we introduce a 
smooth penalty function to the free energy when the order parameters or compositions 
are outside of this range: 

r /i/' -yf if / < 0, 

G'penalty=<0 if /G [0,1], (15) 

[/i(l-/)'-Kl-/) if/>l- 
Here / can be either an order parameter or a composition and /i and v are positive 
constants that are set sufficiently large to avoid non-physical values of /. 

A similar situation exists with regard to the entropy. Expression was derived 
for c^p /Ci G [0, 1], it is undefined if < 0. If is larger than Cj, is defined 
(although not physical) and it can therefore be used as is. We extend the entropy into 
the unphysical parameter space as follows: 



dc 



k^/{n,Ci) [l + ln [cP/C, 



if <C,£, 

(16) 
licP > de. 



Here e is a small parameter that defines the range over which the entropy is modified. 
Expression (fT^ ensures that ds^^'^/dc and its first derivative are continuous at 

5.2. Equilibrium grooving angles 

In equilibrium, a grain boundary intersects an interface at an angle 9 that is determined 
by their relative energies (i.e., the Young-Dupree angle) jS]: 

cos^^i = ^, (17) 
271 

where the geometry is defined in figure 7gb and 71 are the grain boundary and 
interface energies, respectively. This equilibrium angle should also be maintained in 
the immediate vicinity of the intersection even far from equilibrium, since the quantity 
of matter that must be transported to establish the appropriate angle is infinitesimal. 
Hence, this condition is often viewed as a boundary condition. No such condition is 
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Figure 4. A late-time image of two NiSi grains on a Si substrate (the grey scale 
indicates the fraction the NiSi phase present). The angles shown were obtained as 
arccos[7GB/(27i)] and arccos[7GB/(27s)]- 



applied here, the angle is an outcome of the simulations. In this section, we investigate 
whether the phase-field model yields the correct angular condition. 

Figure m shows the example of a NiSi film on a Si substrate and with a free surface. 
The ratio of the grain boundary to twice the interface energies is 7gb/(27i) = 0.441, 
corresponding to 9i = 63.9° and the ratio of the grain boundary to twice the surface 
energy is 7gb/(27s) = 0.235, corresponding to 9^ = 76.4°. The angles drawn in figurelH 
correspond to these theoretical angles. Clearly, the phase-field simulations accurately 
reproduce the theoretical predictions of (fT7|). 

5.3. Platinum segregation 

Figure El shows the Pt concentration profiles in two Ni(Pt)Si grains of a thin film on a Si 
substrate at different times. As Fgb 7^ segregation occurs. The platinum composition 
at the interface, three times as large as the nominal composition, is established very 
early in the simulations. A similar approach can be employed to model Pt segregation 
to other interfaces, e.g., choosing a non-zero value for Isi-msi and ^vacuum-NiSi (i-e., Pt 
concentration-dependent Xsi_NiSi and Xvacuum-NiSi)- 

Figure El shows that the angles are established early in the simulations and are 
maintained even when the surfaces and interfaces do not have constant curvatures (i.e., 
away from equilibrium) as proposed in section 15.21 

6. Conclusion 

We presented a phase-field model for grain boundary grooving in polycrystalline thin 
films on a substrate. The model is quite general, describing compounds and solid 
solution alloys with sufficient freedom to choose solubilities, grain boundary and 
interface energies, and heats of segregation to all interfaces. We also derived analytical 
expressions for describing interface profiles — with and without segregation — and we 
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Figure 5. The Pt composition field at t = 10, 100, 500, and 2000 for the case of 
Ygb = 43 (i.e., Xq-q = 12 — 43cpt). The other parameters used for both simulations 
were Asi = 10, \m = 3, Apt = 3, KNiSi = 2, Kvacuum = 5, kBT/il = 0.4, Xsi-NiSi = 16, 
^NiSi-vacuum = 20 and ^ = 20 for all phases. 



confirmed these using numerical simulation. The present model was shown to accurately 
reproduce the theoretical grain boundary groove angles both at and far from equilibrium. 
As an example, we applied the phase-field model to the special case of a Ni(Pt)Si thin 
film on an initially fiat Si-substrate. The present model is suitable for describing grain 
boundary grooving and thin film agglomeration in real systems. 
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Appendix: grain boundary profile 



In this appendix, we present the derivation of the composition profile at a grain 
boundary. We consider the case of a grain boundary between two semi-infinite grains, 
ijj^^^ and ijj^^^ in one dimension. Since ip^^^ + il)^"^^ — 1 (only two grains of the 
same phase are considered) the only unknowns are the phase fraction of grain 1, 
il>{x) :— ip^^^x) — 1 — ip^'^^x), and the compositions. We assume that there are two 



species on each sublattice i and that cf^^ + cf' = Ci for all i. Hence, all compositions 
are fixed in terms of one composition profile per sublattice, Ci{x). The total energy is 



(2) 



G (V', {C^}) 



dV, (A.l) 



where Yi controls the segregation to the grain boundary for atoms of type i. 
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The functional derivatives of the free energy with respect to if) and Cj are 

^ = 2 I 2D + Xgb - > >^ Q U(l - V^) (1 - 2V^) - 4 «: ^" (A.2) 



In 77^^ Y^i^il-i^-2(>^ + )^)4. (A.3) 



The only dependence of 5G/5ip on the compositions is from Yi Ci. If the deviation of 
Ci from the far field composition is small, then Yi Ci{x) ~ J2i ^-i Ci(oo) and 

^(a;)~l(^l + tanh^), (A.4) 

with 

If there is no segregation, all Y^ are zero and we recover ()12|) and ()13p. 

The equilibrium on sublattice z, given by bGjbci = 0, is independent of the other 
sublattices. We can therefore limit ourselves to one sublattice without loss of generality. 
We will call c the composition of species 1 on this sublattice. If the deviation of this 
composition profile from the nominal composition cq is small, we can describe the profile 
as c(x) = Co + 7(3;) with |7(a;)| << cq. Uding the expression from ()A.4j) for ip in flA.3|l . 
to first order in 7, 6G/6c = gives 

2 (,<.. , _ « -^oosh-^. (A.6) 

This equation is of the form W'^'y"{x) — cP'^ix) = — 5 cosh~^(x/iy), where a and 5 are 
dimensionless constants. The solution to this equation is 

^ _2 ^ ' - aa; —ax 



c{x) ~ Co + - cosh — + ABexp — + A{1 - B) exp + 



r\ +00 9 . 

2rY^/ " -4 

-5y i-irk- 



3 ' a'^-Ak^ 

k=l 



„ 2kx ,^ —2kx 
B exp h 1 — -D exp 



(A.7) 



where A and -B are integration constants. 

Since the general solution given in ()A.7|1 is complicated we rewrite it for two 
particular values of a. For a = 1 

1 ,oX TT ,X ^ , X \ ( , X 



C{X) ~ Co + - 

and for a = 2 



(A.8) 



5 , _2 X 



c(x) ^ Co + - cosh . (A.9) 
W 

Equation p4|) is written for a = 2. Figure El shows il){x) from p2|) and c(x) from (fT^ . 
i.e. flA.9|) . as functions of a;/!^. The circles correspond to simulations using a set of 
parameters such that a = 2. 
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